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Abstract. When the vacuum Einstein equations are cast in the form of hamihonian 
evolution equations, the initial data lie in the cotangent bundle of the manifold TWE 
of riemannian metrics on a Cauchy hypersurface E. As in every lagrangian field 
theory with symmetries, the initial data must satisfy constraints. But, unlike those 
of gauge theories, the constraints of general relativity do not arise as momenta of any 
hamiltonian group action. In this paper, we show that the bracket relations among 
the constraints of general relativity are identical to the bracket relations in the Lie 
algebroid of a groupoid consisting of diffeomorphisms between space- like hypersurfaces 
in spacetimes. 

A direct connection is still missing between the constraints themselves, whose def- 
inition is closely related to the Einstein equations, and our groupoid, in which the 
Einstein equations play no role at all. We discuss some of the difficulties involved in 
making such a connection. 



1. Introduction 

The vacuum Einstein equations state that the Ricci curvature Ric((7) of a lorentzian 
metric g is identically zero. Recast as evolution equations, they become a hamiltonian 
system on the cotangent bundle of the manifolcQ AiT, of smooth riemannian metrics 
on a manifold S which represents the typical Cauchy hypersurfac^. Each element of 
T*AiT, may be identified with a pair (7, tt), where 7 is a riemannian metric and tt is a 
symmetric covariant 2-tensor fieldJl 



2010 Mathematics Subject Classification. 83C05 (Primary); 58II05 (Secondary). 
Key words and phrases. Einstein equations, constraints, Poisson brackets, groupoid. Lie algebroid, 
diffeology. 

Partial support for Christian Blohmann was provided by a Marie Curie Fellowship from the Euro- 
pean Union, for Marco Fernandes by a CAPES grant from the Government of Brazil, and for Alan 
Weinstein by grants DMS-0204100 and DMS-0707137 from the US National Science Foundation. 

""^In this paper, we will actually treat spaces of smooth functions as diffeological spaces rather than 
as Frechet manifolds. In Appendix A, we review the theory of diffeology, concentrating on the aspects 
which are relevant to our work. 

^In general relativity, E has dimension 3, but that assumption is not necessary for anything we do 
in this paper. 

^ Strictly speaking, a cotangent vector tt to A^E should be a contravariant symmetric 2-tensor 
density, but we may use the metric 7 and its associated volume element to identify covariant and 
contravariant tensors, and to identify scalar functions with densities. In addition, we consider tt as an 
endomorphism of the tangent bundle TE when we form tt^ and take traces. 
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It has long been knowr0 that, for a given 7 and n to be admissible as initial con- 
ditions for the Einstein equations, they must satisfy a system of constraint equations. 
These equations may be derived either geometrically from the Gauss-Codazzi equations 
relating the intrinsic and extrinsic curvatures of a hypersurface, or from a lagrangian 
formulation of the Einstein equations in terms of the Einstein-Hilbert action. The 
equations are: 

(1) Cmo(7,7r) := -2div^7r = 0, 

(2) Cen(7, vr) := -R{l) + Ti,{n') - ^^^^(Tr^ = • 

The momentum constraint, C^o, maps T*A^S to the space XT, of vector fields on S, 
while the energy constraint, Cen, takes values in the space J-'S of scalar functions on 
E. The constraints may be viewed as those components of the Einstein tensor which 
involve directions normal to the Cauchy hypersurface. (See Appendix E of ^46] for 
details.) 

The constraint set C C T*Ai'E, where the constraints are all equal to zero, has two 
properties which always hold for the zero sets of momentum^ maps of proper hamilton- 
ian group actions. 

• The constraint set is coisotropic; i.e., for any two functions vanishing on C, 
their Poisson bracket vanishes there as well. (This follows from the bracket 
formulas (jl]) below.) Consequently, on its regular part, namely those pairs 
(7, vr) with no common infinitesimal symmetries, C is foliated by characteristic 
submanifolds whose tangent spaces are spanned by the hamiltonian vector fields 
whose hamiltonians vanish on C. In the group action setting, these would be 
the orbits of the symmetry group. 

• As shown by Arms, Marsden, and Moncrief [1] , the constraint set has quadratic 
singularities at the points which do admit infinitesimal symmetries. In the group 
action setting, this would follow from the linearizability of proper actions and 
the equivariant Darboux theorem. 

The aim of the research described in this paper has been to identify the symmetry 
structure responsible for the constraints and their Poisson bracket relations. What we 
have found is that the bracket relations, rather than coming from the Lie algebra of 
a symmetry group, are those of a Lie algehioid which is derived from a groupo^fQ of 
diffeomorphisms between space-like hypersurfaces in Lorentz manifolds. Unfortunately, 
this groupoid, which encodes the arbitrariness in the choice of initial value hypersurface 
for the Einstein equations as well as of coordinates on this hypersurface, lives over a 
much larger space than the one where the constraints are defined. It remains to be 



''See Section m below for historical remarks and references. 

^Note that we are using two meanings of "momentum" in this discussion, first in a slightly extended 
version of the usual "mass times velocity" , and second in the sense used in the theory of hamiltonian 
actions. In the latter sense, the term "moment" is often used instead. (See the footnote on p. 133 of 
[33] ) for some remarks on the two nomenclatures.) 

^The definition of groupoid is reviewed briefly in Section IA.21 We refer to i^32j for a full treatment 
of Lie algebroids and Lie groupoids. 
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seen what the relevant structure is which connects our groupoid with the constraints 
themselves. 

Since the constraints Cmo and Ccn are vector valued, we break them into scalar- 
valued components to form their Poisson brackets. Following DeWitt [H], we get 
these components by pairing the constraints by integration against vector fields and 
functions on S, obtaining for each vector field X and function (p the following real- 
valued constraint function on T*A^S: 

(3) C(x,<p){l,T^) = / {7(^,Cmo(7,7r)) + 0Con(7,7i")}vol^. 

Since C(x,</,) is the sum of Cx '= C(x,o) and =^ C*(o,(/,), it suffices to write down the 
Poisson bracket relations among these terms. These were found by DeWitt to be: 

{Cx, Cy} = C[XX] 

(4) {Cx,C^} = Cx.^ 

where the metric 7 is here considered as a bundle map from TS to T*S, so that its 
inverse takes the differentials of functions to their gradients. 

The coisotropic property of the constraint set follows immediately from (jl]) above: 
the bracket of any two constraint functions vanishes on the constraint set. On the 
other hand, the dependence of {C^, C^} on the metric 7 means that the brackets are 
not those of a fixed Lie algebra structure on XT, © J-'S. Of course, we may freeze the 
metric 7 to some value 7, but then the resulting bracket { , }y will not satisfy the 
Jacobi identity. The anomaly appears in the jacobiator of a momentum constraint and 
two energy constraints, namely: 

{Cx, {Cc/i, C^}^}^ + circ.perm. = Cjr^(^-i-)(^^dip-iijd(i>) , 

which vanishes only in special cases, such as when X is a Killing vector field for 7. 

We must therefore renounce the idea that the constraints should be the momentum 
map of a symmetry group (such as the group of diffeomorphisms of spacetime). Instead, 
we use a groupoid. 

The objects of our groupoid will be the isometry classes of embeddings of S as a 
cooriented, space-like hypersurface in a lorentzian manifold. We call these S-universes. 
The morphisms of the groupoid, which we call S-evolutions, will be the isometry classes 
of pairs of such embeddings into the same target lorentzian manifold; they can be 
identified with diffeomorphisms between the image hyper surf aces. Remarkably, the Lie 
algebroid of this groupoid is naturally a trivial bundle, and the bracket relations among 
its constant sections turn out to reproduce precisely the bracket relations (jll) among 
the constraints. 

To facilitate computations with the rather abstractly defined spaces of isometry 
classes, we will use the fact that each E-universe has near the image of S a unique 
gaussian representative, i.e. a metric near (in the appropriate sense) S x {0} on S x R 
for which the paths t 1— (x, t) for each x G S are time-like geodesies parametrized by 
(negative) arc- length and normal to S x {0}. 
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As an aside, we will show that the groupoid of S-evolutions is equivalent in a precise 
sense to the groupoid of isometries between those Lorentz manifolds which admit a 
cooriented, space-like hypersurface diffeomorphic to S. 

There remains a significant gap between our results and a satisfactory geometric 
explanation of the relations (jl]) among the constraints, since our Lie algebroid lies over 
a space much bigger than the phase space T*A^S on which the constraints live. The 
latter space may be identified, via a natural riemannian metric on A^S, with the tangent 
bundle TAiH, and this tangent bundle may be in turn identified with the space of 1- 
jets around S of S-universes. Differentiation by tangent vectors based along E gives a 
natural projection from the S-universes to the 1-jets. Unfortunately, it does not seem 
to be possible to make our Lie algebroid descend along this projection. 

To end this introduction, we note that Teitelboim [IS] already gave an argument 
leading to the Poisson bracket relations among the constraints using pure geometry, 
without any appeal to Einstein's equations. In some sense, the accomplishment of 
our paper is to put Teitelboim's argument in its proper mathematical context, that of 
groupoids and Lie algebroids. 
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2. Universes 

We turn now to the construction of a groupoid whose Lie algebroid is a trivial bundle 
with fibre XT, © J^S, such that the bracket of constant sections is given by (jlj). 

2.1. S-universes. A connected lorentzian manifold M will be called a spacetim^^ 
(whether or not its dimension is 4). M will be called S-adapted if it admits a cooriented 
(i.e. with an orientation of the normal bundle), proper embedding of S as a space-like 
hypersurface; such an embedding will be called a S-space in M, and a pair consisting 
of a spacetime and a S-space in it will be called a S-spacetime. 



'Note that any such manifold satisfies an Einstein equation of the form Kic{g) = T if the energy- 
momentum tensor T is simply defined by that equation. 
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Since the identity of specific points in the ambient spacetime is irrelevantly to the 
evolution of metrics on a manifold S, it is natural to make Definition 12.21 below. 

Notation 2.1. For better compatibility with composition, we will often represent map- 
pings by arrows going from right to left; hence we write C°°{X,Y) for the smooth 
mappings to X from Y. 

Definition 2.2. A S-universe is an equivalence class of T,-spacetimes, where M i^H 

i' if) 

and M' ^ S are equivalent if there is an isometry M' M which is consistent with the 
coorientations and which satisfies ip oi = i' . We will denote the se^ of all ll-universes 
by WS. 

To define a diffeology on WE, since different S-universes may be represented by differ- 
ent lorentzian manifolds, we follow the pattern described in Appendix A for functional 
diffeologies on spaces of mappings with varying domains. Namely, we stipulate that a 

parametrization WE P is smooth if: 

(1) each point of P has a neighborhood U for which there is a fixed manifold M such 

that the S-universes p G U, are represented by S-spacetimes Mp -f-> S, 

where all the Mp are open subsets of M; 

(2) ip^p is a smooth map C"^(M, S) ^ P. 

(3) Mtot = {{x,p)\p e P and x e Mp} is open in M x P; 

(4) the lorentzian structures on the fibres are the restrictions of a smooth section 
over Mtot of the bundle of symmetric fibrewise 2-forms; 

The following proposition gives a first description of the tangent bundle of WS. A 
more explicit description will follow the introduction of gaussian splittings. 

Proposition 2.3. The tangent cone to UT, at a T,-universe [{M,g) T] is a vector 
space which may be identified with the quotient of the space r{S'^{T*M)) of symmetric 

2-forms on M by the image of the map T{S'^{T*M)) A'j(s)M taking each vector field 
Z defined on M and vanishing on to the Lie derivative Czg- 

Proof. A tangent vector to is represented by a 1-parameter family {Ms,gs) ^ S of 
embeddings, where Mq = M, Mg are open subsets of a fixed manifold M, and gg are 
lorentzian metrics depending smoothly upon s. 

We will simplify the representative by fixing Mg and the embeddings, so that it is 
only g which varies. To do this, let C,s be the vector field dig/ds = C,g along ig. Since ig 
is a proper embedding for each s, we may "extend" the family ^g over M, i.e. we may 
choose a smooth family ag of vector fields on M such that ^g = ag o ig for all s. This 
family may not be complete, but we may integrate it as far as is possible. The result 

is a family of open sets C M and a family of open embeddings Mg M'^ defining 

^Dirac |16| wrote, "I am inclined to believe . . . that four-dimensional symmetry is not a fundamental 
property of the physical world." Pirani [5^, reviewing Dirac's paper, "finds it difficult to concur". 

^The collection of all S-spacetimes is, like the collection of all sets, a "class" in the set-theoretic 
sense rather than a set. But the S-universes UT, do form a set because every connected manifold is 
diffeomorphic to a submanifold of some M" . 
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a smooth path in C°°(M, [M])open such that Mq = Mq, ipo is the identity on Mq, and 

dips/ds = (Tg o Tpg. Since = o z^, ipg o io = is- (In particular, contains zo(^) 

i' 

for all s.) We obtain a new family (M^, (7^) ^ S of embeddings by setting g'^ = ipl^Qs) 
and i'^ = is o = i^. This family is not quite equivalent to the original one, since 
the domains M'^ are not mapped by ips to all of M^, but they do agree at s = 0. 
According to the observation in Appendix A that tangent vectors are insensitive to 

variations of domain, this insures that the paths in represented by {Ms,gs) ^ S 

and {Mg, g'^) ^ E represent the same tangent vector to WE at [(M, (7) -A S], where the 
latter path has constant embedding i^. We may keep M'^ constant as well, in which case 
the only object varying with s is the lorentzian metric gs- Its derivative with respect 
to s at s = is the required smooth section of S'^{T*M). 

Our representative of the tangent vector is not unique, since we may apply an arbi- 
trary family of diffeomorphisms of M which fix i(E) and the normal orientation there. 
(This freedom corresponds, essentially, to the freedom which we had in extending 
to (Ts.) The only such diffeomorphism which can preserve a lorentzian metric on (con- 
nected) M is the identity. The infinitesimal action of these diffeomorphisms, i.e. the 
map C^g-<r-iC,, is therefore injective, and the tangent space to UT, at [i] is the quotient 
of S^(T*M) by the image of that map. □ 

2.2. S-evoIutions. The relative positions of pairs of E-spaces in the same spacetime, 
up to equivalence, form a groupoid over WE which will be our fundamental symmetry 
structure. 

Definition 2.4. A S-evoIution is an equivalence class of pairs {ii,io) of E-spaces in 
the same spacetime, where a pair {ii, io) in M is equivalent to (i'l, io) if there is a 

single isometry M' i— M which is consistent with the coorientations and which satisfies 
both ip oil = i'l 0''>T'd ip ° io = i'o- denote the set of all E-evolutions by SE. 

The S-evolutions form a diffeological groupoid over the S-universes with the diffeol- 
ogy on SE, like that on UE, defined in terms of representatives. The groupoid structure 

has as target and source the projections [ii,io] [ii] and [ii,io] <— [io] (square brack- 
ets denoting equivalence classes); the composition law is [^2, "^i] [^1, *o] = [^2,'^o]; and the 
inversion rule [21,^0]""^ = [^o,"^!]- 

We will show in Section 12.71 that SE is Morita equivalent to the isometry groupoid 
IE of S-adapted spacetimes. This implies that the orbits of SE are in bijection with 
the orbits of IE, which are the isometry classes of S-adapted spacetimes. 

The isotropy group of [i] consists of all pairs ([21,^0]) such that [ii] = [io] = [i]. For 

such a pair, M E and M E are equivalent, which means that there is an isometry 
ip from M to itself such that ii = ip °io- Such an isometry is unique, if it exists, which 
implies that the isotropy group of [i] is isomorphic to the isometry group of its target 
spacetime. 

Remark 2.5. Each {ii,io) corresponds to a diffeomorphism ii o i^-^ between E-spaces 
in the same spacetime, and composition in SE corresponds to composition of diffeo- 
morphisms. 
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Elements of the Lie algebroid of EYj are infinitesimal S-evolutions and may be 
parametrized by triples consisting of S-universes, "shift" vector fields on S, and "lapse" 
functions on S, per the following proposition. 

Proposition 2.6. The Lie algebroid AST, WE is isomorphic as a vector bundle to 
the trivial bundle WE x (XE © J'E). 

Proof. An element of A£E at the base point [io] G UE is a tangent vector to a smooth 
path in SE whose image is contained in the r-fibre of [io] starting at the unit [io,io]- 
Such a path is represented by a smooth family {ig, i^) defined on an interval containing 
s = 0, and its tangent vector at s = corresponds to a vector field along io. (Note 
that the ig may all be chosen to have the same target; the Lie algebroid fibres are thus 
simpler to analyze than the tangent spaces to UE to which the anchor projects them.) 
Using the lorentzian metric on the target manifold to decompose this vector field into 
its tangential and normal components, and dividing the latter by the unit future normal 
field n, we obtain a vector field X and smooth function on S which correspond to 
the Lie algebroid element given by the path. 

Each fibre of ASE is now identified with the fixed space XE (B J^E. We omit here 
the verification that this identification depends smoothly on the base point in UE. □ 

2.3. Gaussian normal form. To compute the bracket and anchor of the Lie algebroid 
ASE — i- UE, we may work in a neighborhood of the units [i, i]. There, we may use the 
simplified representation of each S-universe [i] near S = i{E) given by the following 
gaussian normal form. Using the metric g on the target of the cooriented embedding 

M <f— > S, we first extend the unit future normal field along 5* by parallel translation 
along the geodesies normal to S to obtain a time-like vector field n on a neighborhood 
U of in M; this extension will satisfy the equations g{n,TS) = 0, ^'(n, n) = —1, and 
Vnii = 0, where V is the Levi-Civita connection of g. Transporting S along the fiow 
of n produces a codimension-1 foliation on a neighborhood of S with space-like leaves 
which are everywhere orthogonal to n. We call the leaves of this foliation the gaussian 
time slices. This construction induces a canonical isometry to a neighborhood of S in 
M from a neighborhood of S x {0} in S x M on which the metric has the gaussian form 

(5) -(7y(x, t)dx^dx^ — dt^) , 

where are coordinates on E. Replacing the former neighborhood by the latter, we 
have established the following normal form result. 

Proposition 2.7. Every E-universe has a representative M -f— > S m which a neighbor- 
hood U of i{E) in M is equal to a neighborhood ofEx {0} in E xW, i{x) = {x, 0), and 
the metric on U has the gaussian form ([5]). This gaussian metric is uniquely determined 
by 

Similarly, any tangent vector to UE has a unique representation in the form 

(6) -aij{x,t)dx^dx^ . 
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When S is compact (and sometimes when it is not), one can take f/ to be a prod- 
uct Ex/ for some open interval / containing zero. We will call such a E-universe 
cylindrical. The "spatial" component ^'yij{x,t)dx'^dx^ of the gaussian metric is then a 
1-parameter family of riemannian metrics on S. Thus we may parametrize the cylin- 
drical S-universes by paths of metrics, and the tangent vectors to them by paths of 
symmetric covariant 2-tensors (not necessarily positive definite) on E. 

For noncompact S, we may have to take U to consist of pairs {x,t) with \t\ < e{x) 
for a smooth positive function e on S. In any case, on a neighborhood of any compact 
subset of S, a S-universe is defined by a path of metrics and a tangent vector to the 
S-universes by a path of symmetric covariant 2-tensors. 

2.4. Gaussian vector fields. To use the gaussian normal form in our computations, 
we must deal with the fact that a slicing which is gaussian for one embedding of S is 
not gaussian for most others. The following lemma about diffeomorphisms will lead us 
to its infinitesimal version about vector fields, which is all we need to use. 

Lemma 2.8. Every diffeomorphism 5 — ?■ S" between space-like, cooriented hypersur- 
faces in spacetimes {M,g) and {M',g') extends to a diffeomorphism ip : U ^ U' between 
neighborhoods of S and S' respectively which respects the gaussian time- splittings, i.e. 
which intertwines the (local) gaussian time flows: 

(7) ^o$f = <o^. 

The diffeomorphism is unique up to the choice of (connected relative to S) U. 

We note that ([7]) holds if and only if ip preserves inner products with the unit normal, 
i.e. g{n,w) = (■?/'* (7') (n, u;) for all vector fields w. By letting {M,g) = {M',g') and 
defining ip in the domain of a flow <l>^ generated by some vector field v on M, and 
differentiating with respect to s, we obtain the following infinitesimal version of Eq. ([7j). 

Definition 2.9. Let U be a neighborhood of a hypersurface S as in Lemma \2.8[ A 
vector field v on U is called ^f-gaussian if it satisfies 

(8) {C,g)in,w) = Q 
for all vector fields w. 

The following infinitesimal version of Lemma 12.81 will be proven by a purely infini- 
tesimal computation. 

Proposition 2.10. Every vector field vq with values in TM defined on a hypersurface 
S as in Lemma \2. <^ extends to a g-gaussian vector field v defined on a neighborhood of 
S. 

Proof. Condition ([H]) can be rewritten as 

= ir,C^g = {C^in + i[n,v\)g = {div + ivd)ir,g + i[n,v\g 
= d{iJng) + iyr,,v]g , 
where we have used the fact that i^g = —dt so that di-^g = 0. Equivalently, 

(9) [n,v] = -gmdg{g{n,v)) , 
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where grad^ is the gradient with respect to g. 

Sphtting the vector field v = X + 0n into components orthogonal and parallel to n 
and observing that [n, v] = VnV — Vi,n = Vnf — Vxn because Vnn = 0, and that 
g{n,v) = —(f) because g{n,n) = —1, we obtain the equivalent condition 

VnW = Vxn + gradg0. 

If we now take the inner product with n, use the facts that gCVnX, n) = and 
5f(Vn0n, n) = — n • (p, then we obtain the condition n ■ = 0. In in other words, 
the gradient grad^0 does not have a normal component, i.e. grad^0 = grad^0, where 
the "spatial gradient" grad^ (j) is defined as grad^ + (n ■ 0)n. (Note the sign again.) 
This implies that that Eq. (Q splits into components orthogonal and parallel to n as 

(10) [n, X] = grad^ , n-0 = O. 

In local coordinates these equations read 

ax 



^11) — =grad^' 



Using the boundary conditions Xt=Q = Xq and (f)t=o = 4>o, we see that the (yf-gaussian 
extension v = X + 0n of the vector field vq = Xq + ^oH exists and is uniquely determined 
in a very simple way by the initial values of X and 0. □ 

We will denote the (/-gaussian extension of Xq + 0on by Gg{Xo, (po). Furthermore, 
we will abuse notation by omitting the zero subscripts for the initial values when the 
context distinguishes them from their extensions. 

For future reference, we write below an explicit formula for the gaussian extension, 
to first order in t. 



(13) GgiXo, 0o) = ^0 + t grad^o '^'o + 4>on + 0{t^). 

2.5. The action of (jf-gaussian vector fields on symmetric 2-forms. To compute 
the anchor ASJ] — )■ TUH of our Lie algebroid, we need to express, in terms of the 
space/time splitting on the ambient manifold, the Lie derivative of a symmetric 2- form 
of the type 

a = -aij{x,t)dx'''dx-^ 

by a (yf-gaussian vector field v = X + 0n. 

The pull-back of the Lie derivative on the ambient manifold to the gaussian time 
shces is 

£ja := C^a — {inCya)dt = C^a — {i[n.v]Cii)dt . 

If V = X + 0n is (^-gaussian we have [n, v] = grad^ 0. Writing = 'C.^ct = 

\(t) Q.ijdx'^dx^ =: d, we obtain 

C^a = Cl^a + 0d + (igrad 4,(y)dt , 
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where Cj^a = £j^(^^a(t) is the Lie derivative on the time shce at t. We will drop the 

superscript of Cj. if it is clear from the context what Cx denotes. 

Using the last equation, we can compute the action of v on the metric, 

CyQ = £^(7 - ^df) = Cxi + 07 + (Vad^07)c?^ - dcpdt 
= Cxi + 07- 

Note that the terms containing derivatives of cancel. Furthermore, we have 

Cng = -2K, 

where K{X, Y) = g(VxY, n) is the second fundamental form with respect to n of the 
time slice at t. (See Section 9.3 of [46], but note that Wald's "extrinsic curvature" is 
the negative of the second fundamental form.) Thus, the Lie derivative of the metric 
with respect to a (7-gaussian vector field is given by 



2.6. The bracket and the anchor. We can now make explicit the Lie algebroid 
structure on ASTj — ?• WE in terms of the trivialization ASH ~ WE x [XE © J^E) given 
by gaussian extension. 

Using (ITUi) . we find the Lie bracket of two ^f-gaussian vector fields to be: 

(15) [X + 0n, r + Vn] = [X, Y] + (p grad^ - grad^ + (X ■ ^/^ - F ■ 0)n , 

which corresponds exactly to (jl]). 

As for any bracket of vector fields, the jacobiator [u, [v,w]] + [v, [w,u]] + [w, [u,v]] 
of three (^-gaussian vector fields u, v, and w vanishes. However, the bracket of two 
(yf-gaussian vector fields is in general not (^-gaussian, as the following proposition shows. 

Proposition 2.11. The bracket l[T5\) of two g-gaussian vector fields is not always g- 
gaussian. 

Proof. We have 

'^■n.C[D,w]9 ■£[?), tu]^n5' ~l~ '^[n,[v,w]\9 



In the second step, we have used that i^g = —dt. On the hypersurface E x {0}, we may 
choose X, 0, y, and ip arbitrarily. For X = F = there, the right hand side of the last 
equation is the second fundamental form contracted with ip grad — grad -0, which 



(14) 



Cx+<t>r,g = Cxi + 07 = Cxi - 20-?^ • 



c?(c/([f,w],n)) + i 

-d{X ■ - r ■ 0) + i[^,grad^ ^\g - «[«;,grad^ <i>\9 

d(X ■ 1p Y ■ 0) + (/^tjigrad^ V ^gra.d^ ipCy)g (-^ui^grad^ gmd^ (f)Cw) g 

-d{X -tp -Y ■ (p) + Cydip - C^d(p + Zgrad^ ^C^g - ig^ad^ ^C^g 

hra^d^ d^Yl - 2lpK) - igra.d^,p{Cxl " 2(pK) 
%rad^ (fyCYl ^grad^ tpCxl ~\~ 2i(0 grad^ ip—ip grad^ <j))K . 



is generally not zero. 



□ 
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The anchor A8^ A TWS of our Lie algebroid is given, up to a sign by the action 
computed in the previous section: 

(16) p{X,(f),g) ■= -Cx{t)+<t>ng = -Ccgix,^)! - , 

where g = 7(t) — ^dt"^ is a gaussian metric on S x M and Gg{X,(j)) = X(t) + 0n 
is the (yf-gaussian extension of X + 0n. This Lie derivative represents the change in 
the "appearance" of the metric on the ambient manifold as the "viewpoint" changes 
according to a vector field X + 0n along a space-like embedding. Note that the anchors 
of two constant sections are applied consecutively to a metric as p{X,(f))p{Y,tlj){g) = 
^Gg{Y,iij)^Gg{x,(t>)g, that is, by applying the Lie derivatives of the gaussian extensions in 
reverse order. This takes care of the negative signs in Eq. f|T6|) . 

The kernel of p consists of those {X, 0) whose (7-gaussian extension is a Killing vector 
field. It follows that the dimension of this kernel is zer over an open dense subset of 
WE, since any S-universe is contained in a 1-parameter family whose generic members 
have no isometrics. (Local perturbations of the lorentzian metric suffice to achieve 
this; deeper results in the riemannian case go back at least as far as pJZj-) It follows 
immediately that the bracket on sections of the Lie algebroid ASH is determined by that 
on their images under the anchor (fT6|) . In particular, the bracket of constant sections 
{X, 0) and {Y, ip) of the trivial bundle must be that given by given by Eq. (fT5|) : 

(17) [(X, 0), (y, V^)] = ([X, y] + grad^ ij-ij grad^ 0, X ■ - F ■ 0) , 

where we now evaluate X and at t = 0. Together with the anchor flTBl) . this bracket 
determines the Lie algebroid structure on AUTj. To summarize, we have: 

Theorem 2.12. The T,- evolutions £S form a groupoid over the Tj-universes UT,. Each 
orbit of this groupoid consists of all ll-universes which are represented by T,-spaces in 
a fixed spacetime {M,g). The Lie algebroid AST, has a natural identification with the 
trivial bundle UTx {XT(BJ^T). Under this identification, in the gaussian representation 
of ambient metrics, the anchor is the Lie derivative by gaussian extensions: p{X, 0, g) = 
CGg{x,<j))g- The bracket of constant sections is given by ( (77| ). 

2.7. An equivalent groupoid and the moduli stack of spacetimes. This section 
is peripheral to the main argument of this paper, but it suggests another point of view 
toward the groupoid ST of E-evolutions. 

We have just seen that the orbits of ST are in one-to-one correspondence with isom- 
etry classes of E-adapted spacetimes, while the isotropy groups in ST are just the 
isometry groups of those spacetimes. It turns out that ST is equivalent to another 
groupoid in which T plays a much less central role. We only require that the space- 
times be S-adapted, without ever specifying the placement of T. For the notion of 
equivalence of groupoids, we refer to [TOj and [34j . 

Definition 2.13. The groupoid XT is defined to be that in which the objects are T- 
adapted spacetimes, and the morphisms are isometrics between these spacetimes. 



The dimension of the kernel of p is always finite; it is at most ^(n + l)n + n + 1, where n — dim S, 
with equality only when g has constant sectional curvature. 
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Remark 2.14. Strictly speaking, XS is not a groupoid, since the spacetimes do not 
form a set. But it is equivalent to its wide subgroupoid consisting of isometries between 
those spacetimes whose underlying manifolds are submanifolds of euclidean spaces. 

IE may be identified with an action groupoid whose objects are the elements of the 
collection iSS of all S-adapted spacetimes (or, alternatively, a set as in the remark 
above). The groupoid acting on SH consists of the diffeomorphisms between these 
objects. If {Mi,gi) is a S-adapted spacetime and $ is a diffeomorphism to a manifold 
Mq from Ml, then the result of acting on {Mi,gi) by $ is defined to be the spacetime 
(Mo, ($~^)*5'i). $ then becomes an isometry and so may be considered as a morphism 
in XS. 

Proposition 2.15. The groupoids XT, and ST are equivalent. 

Proof. An equivalence between groupoids is given by a biprincipal bundle, i.e. a space 
on which the groupoids have commuting free actions, with the fibres of the moment 
map0 of each one being the orbits of the other. Such a bibundle induces bijections 
between the orbit spaces of the two groupoids and between the isotropy groups of 
corresponding objects. 

To get an equivalence between XS and ST, we take as total space of our bundle the 
collection "HE of all S-spaces in all possible spacetimes. 

The moment map for the left action of XT forgets the embedding and remembers only 

the target. The typical morphism in XT is an isometry M' ^ M between spacetimes. 

i ipoi 

It acts on any S-space M <(— > S to give the (equivalent) embedding M' M. This 
action is free because an isometry of a connected manifold which fixes a hypersurface 
and its normal bundle must be the identity. 

The moment map for the right action of ST takes each M ^ M to its equivalence 
class [i]. If i G TiC and [ii,io] are such that [i] = [ii], then we define the right action 

i i' 

of [ii, io] on i as follows. The equivalence between M > S and M' ^ S is realized by 

a unique isometry M' M such that ip o ii = i. Then i ■ [ii,io] is defined to be the 

embedding M ^ T. To see that this action is free, suppose that i ■ [ii,io] = i- Then 
M = M' and ip o = i = ^ip o i^, so ii = io, and [ii, io] is an identity morphism. 

The transitivity of the left action on the right moment fibres is just a restatement of 
the definition of equivalence used in defining the E-universes. Transitivity of the right 
action on the left moment fibres is obvious, since the morphisms in ST are (equivalence 
classes) of pairs of embeddings into the same target. □ 

Remark 2.16. The notion of stack was introduced in algebraic geometry and has 
recently migrated to differential geometry [6]; the purpose of the notion is to provide 
a description of spaces of equivalence classes when it is important to keep track of the 
multiple ways in which objects can be equivalent and thereby to overcome difficulties 
related to singular behavior of quotient spaces. One way to understand stacks is to see 

-'^^An action of a groupoid G over Gq on a space X includes as part of its data a map Gq -i^ X 
which determines which groupoid elements act on which elements of X. This map is sometimes called 
the moment map of the action. 
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them as presented by groupoids, where equivalent groupoids determine "the same" 
stack in the same sense that a given manifold may be described in different ways by 
overlapping families of coordinate charts. 

The equivalence between ST, and XS which we have just proven shows that we may 
consider ST as a presentation of the moduli stack of S-adapted spacetimes, i.e. the 
isometry classes (including information about self-isometries) of spacetimes admitting 
S as an "instantaneous space", or "initial condition". (Note that this is a purely 
"kinematic" construction, as we have not imposed any dynamical condition such as the 
vacuum Einstein equations.) 

3. Discussion: the descent problem 

We have constructed a groupoid over UT whose Lie algebroid bracket, for constant 
sections in a natural local trivialization, exactly matches the bracket relations on the 
constraints for Einstein's equations. To establish a more direct relation with the con- 
straints themselves, we would need to find similar structures on the phase space T*AiT. 

There is a natural projection V to TAiT from UT, assigning to every S-universe the 
1-jet with respect to t at t = of its gaussian representation as a path in AiT. (Note 
that this is well defined even if the S-universe is not cylindrical.) But there is no way 
to push our Lie algebroid forward under this projection, essentially because the value 
of the anchor at a given 1-jet would have to depend on the 2-jet. 

To surmount this difficulty, we tried to use a second-order evolution equation on 
AiT, i.e. a rule which expresses 2-jets in terms of 1-jets, such as the Einstein evolution 
equations themselves. But the resulting anchor was not consistent with the bracket 
relations. 

We also tried "reverse engineering" , defining the anchor so that it would take constant 
sections to the hamiltonian vector fields of the constraint functions. But the anchor 
so-defined, when applied to two constant sections associated to functions and on 
M, generally takes their bracket as defined by the metric-dependent relations (jl]) to a 
vector field which is not even hamiltonian. 

4. Some history 

A hamiltonian formulation of general relativity can be found in the work of Pirani and 
Schild [37] , aimed at the quantization of Einstein's gravitational field equations. These 
authors, as well as Bergmann, Penfield, Schiller, and Zatzkis [9] consider a physical 
state at a certain time to be given by data on a space-like hypersurface, which must in 
some sense be arbitrary in order to maintain four- dimensional covariance. 

Bergmann [7| had by then already begun a systematic study of covariant field theories 
of general type, addressing the problem of bringing general relativity into the canonical 
form as a preliminary step to quantization. While canonical quantization of field theo- 
ries was being developed, it soon became clear that general relativity posed additional 
difficulties connected with the degeneracy of the lagrangian which was a consequence 
of four- dimensional diffeomorphism invariance. In fact, Dirac's original work [15] on 
constrained dynamics was inspired in large part by this problem. 
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It was noted by Pirani, Schild, and Skinner [38] and not long after by Dirac fT6] and 
Arnowitt, Deser, and Misner (ADM) [3J that great simplifications could be made at the 
expense of giving up four-dimensional symmetry. Such simplifications became possible 
by fixing a foliation of spacetime by space-like hypersurfaces, with the physical states 
living on these surfaces. Such a spacetime decomposition leads to the decomposition of 
vectors along each hypersurface into their normal and tangential components, and the 
metric tensor itself may be presented in the form 

„(4) + N,N-^ Ar„ \ 

\ -''m ymn I 

where the lapse function and the shift 3- vector A^" describe the variation of the time 
and space coordinates on infinitesimally close space-like hypersurfaces. The lapse and 
shift can be chosen arbitrarily but enter in the constraint equations, which arise from 
the degeneracy of the lagrangian and are given by time components of the Einstein field 
equations G^^ = 0. It had already been realized by Dirac and many others that the shift 
functions A„ generated coordinate transformations on the hypersurfaces, whereas the 
lapse A^ was related to time translation. In this sense, lapse and shift could be viewed 
as gauge potentials that had to be fixed in order to solve the initial value problem. 
The lapse and shift functions were first introduced in [2]. Their geometrical meaning 
was explained in |17] and in more global terms by Fischer and Marsden [18], who were 
perhaps the first authors to suggest that the constraints should be seen as something 
like a momentum map. 

Still in search of avoiding the shortcomings of coordinate-dependent language, Kuchaf 
[25] argued that field dynamics does not take place in spacetime, or along a single 
foliation of hypersurfaces but in what he called hyperspace, an infinite-dimensional 
manifold consisting of all the spacelike hypersurfaces in a given spacetime. 

After the appearance of [16], Katz [26] found the formulas for the Poisson brackets 
of the constraint functions. In their space-integrated form, the brackets were first 
computed by DeWitt [l^ . While the computation of the brackets was straightforward, 
their geometric interpretation was not satisfactory. A number of authors have tried to 
give a more conceptual derivation, e.g., by studying hypersurface deformations [U], as a 
method to guarantee the path-independence of geometrodynamical evolution [27], or by 
generalizing the concept of transformation groups [8] [17] . In particular, Teitelboim [15] 
was perhaps the first to show that the Poisson bracket relations are purely a consequence 
of the geometry of hypersurfaces in a riemannian or lorentzian manifold, independent 
of any particular field theory. In a sense, our paper may be seen as setting Teitelboim's 
argument in its proper mathematical setting, that of groupoids and Lie algebroids. 

The most ambitious approaches aimed at recovering diffeomorphism covariance of the 
initial value problem [21] [21] [25] [2S]- Parallel developments were also made by Wheeler 
[17]. Inspired by the dynamical description of the electromagnetic field in terms of the 
vector potential, he proposed to describe the dynamics of three-space metrics through 
the propagation of the intrinsic metric of space-like hypersurfaces with respect to a time 
coordinate. Such a time parameter would label the leaves the dynamically produced 
spacetime. In this approach the extrinsic curvature of these hypersurfaces corresponds 
to the canonical momenta. The configuration space of this theory is known as Wheeler's 
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superspace. Trajectories of metrics in this space produce four-dimensional spacetime 
geometries. 

Since the bracket of constraint functions is metric dependent, this suggested to 
Bergmann and Komar [8] that the associated symmetries could be metric-dependent 
as well. Their paper contains many ideas which are very close to ours, although the 
language is somewhat different. In particular, we would say that their "Q-type trans- 
formations" are precisely the bisections of the action groupoid (essentially our XS, but 
without a specific choice of S), associated to the diffeomorphisms of a 4-manifold M 
acting on the function space of lorentzian metrics on M. Their infinitesimal transfor- 
mations are the sections of the action Lie algebroid; their equation (3.1) is precisely 
the formula for the bracket of sections in this Lie algebroid! Bergmann and Komar also 
observe that the orbits are isometry classes and that the action fails to be faithful in 
the presence of isometries. 

Bergmann and Komar even make the tantalizing statement, "That these transforma- 
tions form a group, or at least a groupoid, is seen from their definition." Unfortunately, 
groupoids do not reappear anywhere in the paper, and it is not clear what notion of 
groupoid the authors had in mind. In particular, although [8] includes a discussion of 
the set of diffeomorphisms between hypersurfaces in spacetime (the infinitesimal trans- 
formations being 4- vector fields along these hypersurfaces), there is no suggestion that 
these form a groupoid. 

The idea of associating the Q-type transformations with diffeomorphism invariance 
of general relativity appeared also in |Tl] where these transformations arise from "field 
dependent" gauge generators. These field dependent generators appeared also in [39] . 
110], [S], and [12], and more recently in |35j . 

Hojman, Kuchaf, and Teitelboim |2T] look at space-like embeddings into a fixed 
lorentzian manifold and find what is more or less the Lie algebroid bracket for the Lie 
algebroid of the groupoid of diffeomorphisms between hypersurfaces. 

Finally, we should at least mention the immense analysis literature on existence 
theory, both local and global, for the Einstein initial value problem, beginning with 
fundamental work of Lichnerowicz [31j and Foures-Bruhat (= Choquet-Bruhat) [IH] 
and continuing to this day. (See [13] for a fairly recent survey.) We hope that our work 
provides new geometric understanding which may contribute to both the analysis and 
the quantization of the Einstein field equations. 

Appendix A. Diffeology 

Although some of the infinite-dimensional spaces of smooth mappings in this paper 
may be considered as Frechet manifolds, this analytical structure is not necessary for 
formal computations. Instead, we work in the framework of diffeological spaces. These 
objects were introducecQ by Souriau [43j and developed extensively by Iglesias |22] and 
others. We give a brief introduction to diffeology here and refer to [SU], [50], and the 
work-in-progress [23] for further details. Similar notions in the topological setting are 
discussed in [1]. 



-'^^A very similar notion was introduced by Chen [T^j, and different notions of "smootheology" are 
compared in [33] 
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Roughly speaking, we can do a great deal of differential geometry on a set X once we 
know what it means for a family of elements of X to depend smoothly on parameters. 

Definition A.l. A parameter space is an open subset of MJ^ for some n, and a 
parametrization of a set X is a map to X from some parameter space P. A diffeol- 

ogy on X is a set T> of parametrizations of X which contains all constant maps, which 
is closed under composition on the right with smooth (in the usual sense) maps between 
parameter spaces, and which is locally defined in the sense that a parametrization is in 
T) if and only if its restrictions to all the sets in some open covering are. The elements 
ofD are called 'plots , or smooth parametrizations. (X, P) is called a diffeological 
space; we denote it simply by X when it is clear what diffeology is being used. A plot 

X ^ P with X = 0(0) is called a plot at x. 

Remark A. 2. In the language of sheaf theory, a diffeological space is a concrete sheaf 
on the site of open subsets of euclidean spaces [S]. This point of view is particularly 
well suited to show that the category of diffeological spaces has small limits (taken 
point- wise), small colimits (first taken point-wise, then sheafified), and exponential ob- 
jects (given by the universal property), thus allowing for constructions such as quotient 
spaces, pull-backs, mapping spaces, etc. that generally fail to exist for smooth mani- 
folds. For our purposes, however, we will need to give the explicit descriptions of these 
constructions. 

Just as in topology, every set X carries the discrete diffeology, for which only the 
locally constant maps are plots, and the coarse diffeology, for which every parametriza- 
tion is smooth. If X is a (finite-dimensional) manifold, the usual smooth maps to X 
from parameter spaces form the "standard" diffeology. 

Diffeological spaces are the objects of a category in which the morphisms are the 
smooth maps, defined as follows. 

Definition A. 3. A map X between diffeological spaces is a smooth map if f o(j) 
is a plot for X whenever (p is a plot for Y . We denote the set of all such smooth maps 
by C°°{X,Y) . (Note that many authors denote it by C°°{Y,X).) A smooth map with a 
smooth inverse is a diffeomorphism. 

A smooth map between manifolds with the standard diffeologies is just a smooth map 
in the usual sense. The parametrizations of any diffeological space which are smooth 
are just the plots, so the term "smooth parametrization" has an unambiguous meaning. 

Any diffeological space {X, V) carries the P-topology, defined as the finest topology 
for which all plots are continuous; i.e., a subset f/ C A is open [closed] if and only if 
(f)~^{U) is open [closed] for every plot (f). Smooth maps are always D-continuous. 

A product A X y of diffeological spaces carries a product diffeology, whose plots are 
the parametrizations whose compositions with the projections to X and Y are smooth. 
Each subset of a diffeological space has a natural subspace diffeology in which the 
plots are those parametrizations whose composition with the inclusion is smooth. The 
restriction of a smooth map to any subset with the subspace diffeology is again smooth. 

Any set Dq of maps to X from diffeological spaces generates the diffeology T>q 
consisting of all those parametrizations which are locally compositions of the form 0os, 



GROUPOID SYMMETRY AND CONSTRAINTS IN GENERAL RELATIVITY 



17 



where </> is in Vq and s is a parametrization of the domain of 0, together with all 
constant maps from parameter spaces. (The latter are already included if the images 
of the elements of Vq cover X.) 

If consists of a single surjective map X ^ Y from a diffcological space Y, the 
diffeology which it generates is called the quotient diffeology. Conversely, any diffe- 
ology is the quotient diffeology for the union of all its plots, considered as a single map 
defined on the disjoint union of the domains of the plots. 

If X and Y are diffeological spaces, the functional diffeology on C°°{X, Y) is that 

for which a parametrization C°°{X,Y) -f^ P is a plot if and only if the corresponding 
evaluation map X ^ Y x P is smooth. The "exponential law" C°°{C'^{X,Y), Z)) ^ 
C°^(X, Y X Z) then holds for all X, Y, and Z, not just in the defining case where Z is a 
parameter space, and the composition operations C°^(X, Z) C°°(X, Y) x C°°{Y, Z) 
are smooth. 

The subspace diffeology construction produces diffeologies on spaces of mappings 
satisfying extra conditions, such as spaces of diffeomorphisms or embeddings, spaces 
of sections of smooth bundles (e.g. tensors), and solution spaces of ordinary or partial 
differential equations. 

We can also define diffeologies on spaces of mappings with variable domains. For 
simplicity, we let these domains be subsets of a fixed space. Let \Y] be some collection 
of subsets of a diffeological space Y (for instance the open subsets, if Y carries a 
topology), and let C°^(X, \Y]) be the set of all mappings to X whose domains belong 

to \Y\ . Thinking of a parametrization C°° (X, \Y\ ) <(— P as a family X 1^ of maps 
parametrized by p e P, we call it a plot when the evaluation map X-<— DCFxPis 
smooth, where the domain D = {{y,p)\y G 1^)} carries the subspace diffeology. When 
y is a manifold and [Y] consists of the open subsets, it may be appropriate to restrict 
the diffeology to consist of those families for which the domain D is open in F x P, so 
that we are always dealing with maps defined on manifolds. In this case, we denote the 
space of mappings by C°°{X,Yopen)- In particular, if X is a single point, we obtain a 
diffeology on the set of open subsets of Y. 

Example A. 4. If S is any manifold, the space of riemannian metrics on E carries 
a functional diffeology. So does A412open, the metrics defined on open subsets of E. 

If the diffeological space Y carries a topology (not necessarily the P-topology), 
C°°(X, y) also carries the "finer" compact functional diffeology in which the plots 

C°°{X, Y) <^ P are required to satisfy the additional condition that each p E P has a 
neighborhood V for which all the maps in 0(V) agree outside some compact subset of 
Y. C~(X, Y) with this diffeology is denoted by C^{X, Y). 

A.l. Jets and tangent vectors. To define jets in diffeology, we start with the basic 
case of parameter spaces. For k = 0,1,2, ... ,oo, and smooth maps s and t to M" from 
a parameter space P containing 0, s and t are defined to have the same k-jet at if 
s(0) = t(0) and if their partial derivatives through order k match at 0. Plots / and g 
to a diffcological space X from P will be said to have the same A';-jct at if there is a 

plot X Q such that f = ho s and g = hot, where s and t are plots for Q the same 
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/c-jet at 0. Finally, for any diffeological spaces X and Y, two maps X Y have the 
same fc-jet at |/ e F if their compositions with any plot at y have the same /c-jet at 0. 
(It is easy to see that, if X and Y are parameter spaces, this coincides with the original 
definition.) 

Having the same A:-jct at y is an equivalence relation and two maps with the 
same germ at y have the same fc-jet there for any k. We thus obtain an equivalence 
relation on germs, and we define the space J*^(X, F) of /c-jets of maps to X from Y to 
be the set of pairs {jyf, y), where / is the germ of a map to X from some neighborhood 
of y in y, and j^f is the equivalence class of / for 

\i K (ZY , wc say that / ~^ g H f and g have the same fc-jet at all points of K, and 
we call the equivalence classes for this relation /c-jets along k. If X is a bundle over Y, 
we may refer to the k-jets of those maps y ■<— X which happen to be sections of the 
bundle as /c-jets of sections. 

We define diffeological structures on jet bundles by considering them as quotients of 
spaecs of mappings, using the functional diffeology. When k = oo and X and Y are 
manifolds, the sheaf diffeology is also interesting, since the smooth maps are maps into 
the leaves of a fohation. 

Since constant maps are always smooth, and the 0-jets of maps are just their values, 
there is a natural identification of J^{X, Y) with X x Y. There are also natural maps 
J^{X,Y) ^ J''{X,Y) for k < I. For k = 0, this gives natural projections of the jet 
spaces to X and Y. That the infinite jet space J°°(X, Y) is the inverse limit of the jet 
spaces for finite k. 

Jets of mappings into mapping spaces arc mappings into jet bundles. If X, Y, and Z 

are manifolds, smooth maps C°°(X, Y) Z correspond to smooth maps X ■'^Y x Z. 
It follows that J^(C°°(X, Y),Z)= C°°{J^{X, Z),Y). 

Of special importance are the 1-jets at of maps to X from neighborhoods of in 
R; these are the tangent vectors. We denote the set of all tangent vectors to X by TX 
and call it the tangent cone bundle of X. It is the disjoint union of tangent cones 
T^X at the points of X, which are cones because reparamctrization of curves by the 
action of the multiplicative group M\{0} on M leads to a natural action of this group 
on the tangent spaces. TX has a diffeological structure and projection to X inherited 
from those on J^(X,R). 

For mapping spaces between manifolds, the general result above on jet bundles gives, 
with Z = R and k = 1, TC°°{X,Y) = C°°{TX,Y). The tangent bundle projection 
C°°{X, Y) ^ TC°°{X, Y) is just composition with the projection X ^ TX. 

Two warnings are in order. First of all, the action of K\{0} on the tangent cone 
may not be faithful. For instance, if X is the half-line [0, oo) viewed as M/Z2 with the 
quotient diffeology, its tangent cone at is also a half-line on which multiplication by 
— 1 is the identity. Second, it is generally not possible to add 1-jets of curves. If it were, 
we would find that adding a vector v at in R/Z2 to itself would give both 2v and 0. 
Another example where addition is not possible is the tangent cone at to the union 
of the coordinate axes in with the subspacc diffeology. 

An example where all the tangent cones are zero is the space [Fjopen = C'°°({p}, [^])opei 
of open subsets of a manifold Y. A path in [Fjopen is just an open subset of / x F for 
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some interval / of real numbers. If Us is any such path defined around s = 0, it has 
the same tangent vector as the constant path through Uq. In fact, they factor via 
curves tangent at through the plot on / x M defined by f/(s,o) = Us, U(s,s'2) = ^o, 
and U(^s,s') = y elsewhere. Similarly, tangent vectors to spaces C°°(X, [F])opcn are in- 
sensitive to the variation of domain with s and may all be represented by families of 
functions with unchanging domains. The same arguments apply to jets of any order. 

Example A. 5. For the space A^Sopen of riemannian metrics defined on open subsets of 
S, the tangent space at a metric g defined on all of S consists of the smooth symmetric 
covariant 2-forms on E, even though the domain of a path through g may shrink as the 
path parameter varies. 

It is possible to define a tangent vector space at each point of a diffeological space by 
taking formal linear combinations of 1-jets of curves, as in [20J, but it is then necessary 
to introduce further relations, so that, for instance, the tangent space at the conical 
singular point of the space R/Z2 reduces to zero, while the tangent cone does not. 
(On the other hand, the tangent space at the intersection point in the union of the 
coordinate axes in the plane is two-dimensional.) 

For a diffeological group, group multiplication induces a a vector space structure on 
each tangent cone [30] . 

Cotangent spaces to diffeological spaces may be defined as spaces of 1-jets of smooth 
mappings to M; these all have vector space structures derived from that on M. The 
composition of real-valued functions with curves defines a natural pairing (which can 
be degenerate) between tangent and cotangent spaces. 

If N and E are manifolds and p is a submersion, then the finite jet spaces J^{Y, E) 
are also smooth manifolds with submersions to Y , and the infinite jet space J°°{Y,E) 
carries the projective limit diffeology in which a map into it is smooth if all of the 
compositions with projections into finite jet spaces are smooth. 

A. 2. Diffeological groupoids. Diffeological groupoids are defined like topological 
groupoids [22j[23j. Recall that a category C consists of a collection Cq of objects 
and a collection Ci of morphisms, with target and source maps I and r to Cq from 
Ci, a unit inclusion map Ci 4- Cq, and a composition map to Ci from C2 = {{f,g) G 
Ci X Ci|r(/) = 1(g)}, satisfying the usual axioms for associativity and units. A group- 
oid is a categor}o in which every morphism is invertible. 

If Cq and Ci are diffeological spaces, we give C2 C Ci x Ci the subspace diffeology. 
If all of the structure maps are smooth, we say that C is a diffeological category. If C 
is, in addition, a groupoid, and the map l which takes each morphism to its inverse is 
smooth, then C is a diffeological groupoid. Note that we do not require I and r to 
be submersions; in fact, this notion is better replaced with that of "subduction" [23] in 
the diffeological case. But if a diffeological groupoid C is a manifold and / and r are 
submersions, then C is a Lie groupoid in the usual sense. 

-'^"'Many authors require a groupoid to be a small category, in the sense the morphisms and objects 
form sets rather than just cohections. We will not make this assumption, but will occasionally point 
out how to replace "large" groupoids by small ones. 
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Defining the Lie algebroid of a Lie groupoid is not so simple, even wfien Co is a single 
point, in which case Ci is a diffeological group. A "Lie algebra" bracket for such a 
group G is defined in using the conjugation operation of G on itself. It is a bilinear 
operation on the tangent (vector) space TeG at the identity, but antisymmetry and the 
Jacobi identity have been established in [30] under some extra assumptions on T^G, 
holding for example in the case where G is the group of diffeomorphisms of a manifold 
M with the functional diffeology. In this case, T^G is, as expected, the space XM of 
vector fields on M, and the bracket is the usual Lie algebra bracket. 

If we replace M above by a diffeological space X, it is no longer even clear that 
inversion is smooth in the group of diffeomorphisms, though this is the case for many 
X. We can always define a stronger diffeology by admitting as plots only those whose 
composition with inversion is a plot in the functional diffeology. The tangent space at 
the identity may be identified with those vector fields which are tangent to paths of 
diffeomorphisms. 

Let G be a diffeological groupoid, and let B{G) be its group of bisections, i.e. smooth 
sections 7 of r for which / 07 is a diffeomorphism. We give B{G) the diffeology in which 
the plots are those plots for the subspace-functional diffeology for which composition 
with inversion is also a plot. The tangent space to the identity then consists of smooth 
maps TGi A Gq lifting the unit section e which are tangent to smooth paths through 
the identity in B{G). The values of a are tangent to the r-fibres, so it is natural to 
consider them as sections of the "bundle" A{G) over Go which is the puUback by e of 
ker Tr. Without further assumptions, A{G) is not a vector bundle, but it should still 
play the role of the Lie algebroid. To get a bracket operation, we follow [20] and [30] 
and use the natural action of the group B{G) on A{G). Its derivative at the identity 
with respect to the first variable gives a binary operation [a,b] on sections of A{G), 
defined a priori only when a is admissible in the sense of being tangent to a path 
in B{G). In addition, there is an anchor which takes admissible sections of A{G) to 
admissible vector fields, i.e. admissible sections of TGq. 

Since the action of B{G) preserves all the structures in sight, the operation of bracket- 
ing on the left by an admissible section is a derivation with respect to both the bracket it- 
self and multiplication by functions. In other words, we have: [a, fb] = f[a, b] + {p{a)f)b 
for admissible a, functions /, and all sections b, and [a, [b,c]] = [[a,b],c\ + [b, [a,c]] for 
admissible a, b, and c. 
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